Noname manuscript No. 

(will be inserted by the editor) 



On the evolution equations for ideal 
magnetohydrodynamics in curved spacetime 

Daniela Pugliese • Juan A. Valiente 
Kroon 



Received: date / Accepted; date 



(N 



O 

(N 



Abstract We examine the problem of the construction of a first order sym- 
metric hyperbolic evolution system for the Einstein-Maxwell-Euler system. 
Our analysis is based on a 1 + 3 tetrad formalism which makes use of the 
^ \ components of the Weyl tensor as one of the unknowns. In order to ensure 

the symmetric hyperbolicity of the evolution equations implied by the Bianchi 
0^ i identity, we introduce a tensor of rank 3 corresponding to the covariant deriva- 

tive of the Faraday tensor. Our analysis includes the case of a perfect fluid with 
^ I infinite conductivity (ideal magnetohydrodynamics) as a particular subcase. 



in 



Keywords Magnetohydrodynamics, initial value problem 



■ 1 Introduction 



The present article is concerned with the construction of a first order system of 
i/^ ■ quasilinear hyperbolic evolution equations describing a charged ideal perfect 

fluid coupled to the Einstein field equations — the Einstein-Maxwell-Euler sys- 
. tern. As the fluid is charged, one needs to bring into consideration the Maxwell 

equations to describe the electromagnetic field produced by the charged fluid 
flow. This system contains as an important subcase that of the so-called ideal 
magnetohydrodynamics (MHD) — which applies applies to the study of plasma 
in situations where a fluid subjected to significant magnetic fields. 
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As it is well known, General Relativity admits initial value problem for- 
mulation whereby one prescribes certain initial data on a 3-dimensional hy- 
pcrsurface, and one purports to reconstruct the spacctimc associated to this 
initial data — a so-called Cauchy problem. The formulation of an initial value 
problem is a necessary starting point for the construction of numerical solu- 
tions to, say, the Einstein-Maxwell-Euler system. The importance of a initial 
value formulation of the system under consideration is not restricted to nu- 
merical considerations, it is also a natural starting point for a wide variety of 
analytical studies of the qualitative properties of the solutions to the equa- 
tions. In this article we approach the construction of the evolution equations 
of the Einstein-Maxwell-Euler system from the point of view of mathematical 
Relativity. Hence, the amenability of our analysis to analytic considerations 
takes precedence over numerical considerations of implementability. Examples 
of qualitative aspects of the solutions to, say, the Einstein-Maxwell-Euler sys- 
tem requiring an suitable initial value formulation are the discussion of local 
and global existence problems and the analysis of the stability of certain ref- 
erence solutions. Our initial value formulation of MHD will be used elsewhere 
to discuss the stability of certain spherically symmetric configurations using 
the methods of [1]. 

In order to formulate an initial value problem for the Einstein-Maxwell- 
Euler system, one has to break the covariancc of the equations by introducing 
coordinates and choosing a preferred timclikc direction in the spacctimc. This 
choice of coordinates and timelike direction is usually known as a gauge choice. 
In matter models which contain fluids, there is a natural timelike direction — 
that given by the fluid flowlines. The reasons behind the need of a gauge 
choice to formulate an initial value problem are technical: the machinery of 
the theory of partial differential equations does not apply to tensorial objects. 
As a consequence, the discussion of initial value problems in General Relativity 
is inherently gauge dependent. 

The motivation behind the gauge choice procedure discussed in the pre- 
vious paragraph is to deduce, from our basic set of (tensorial) equations, a 
closcqj subsystem of evolution equations which is in what is called symmet- 
ric hyperbolic form — see e.g. [5] for precise definitions. This procedure is 
sometimes called hyperbolic reduction. Once a symmetric hyperbolic system 
of equations has been obtained then the machinery of the theory of hyperbolic 
equations applies and ensures that for suitable data prescribed on an initial 
hypersurface, the evolution equations have a local solution (i.e. a solution that 
exists at least for for a small interval of time) and that the solutions thus 
obtained is unique. The same theory also ensures that the solutions obtained 
from the Cauchy problem depend continuously on the values of initial data. 
An initial value problem with the properties described in the previous para- 
graph is said to be well-posed. It is important to point out that if a certain 
reduction procedure does not lead an hyperbolic system, this does not imply 



^ Closed in this context means that there should be as many equations as unknowns. 
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that the resulting problem is not well posed — this this would require a much 
more detailed analysis. Symmetric hyperbolicity is just a necessary condition 
for well-poscdncss. 

The purpose of the present article is to show that the Einstein-Maxwell- 
Euler system admits a reduction procedure leading to a system of symmetric 
hyperbolic evolution equations. This system of equations is of first order. As 
a consequence of our construction one automatically obtains a local existence 
and uniqueness result for the system under consideration. 

The hyperbolic reduction procedure described in the present article bor- 
rows from the discussion of the evolution equations for the Einstein-Euler 
system by H. Friedrich in [5] — see also [HH] and [S]. In this reference a La- 
grangian gauge was used to undertake the hyperbolic reduction and to obtain 
the desired evolution equations. The term Lagrangian means, in this context, 
that one uses the flow lines of the fluid to construct the preferred time direc- 
tion on which the hyperbolic reduction is based. As it will be seen in the main 
body of the article, this choice has important technical advantages. The central 
equation in this discussion is the Bianchi identity. It provides evolution equa- 
tions for the components of the Weyl tensor. The addition of electromagnetic 
interactions to Friedrich's system destroys, in principle, the symmetric hyper- 
bolic nature of the evolution equation as derivatives of the Faraday tensor 
enter into the principal part of the Bianchi evolution equations. This difficulty 
can be handled by the introduction of a new field unknown corresponding 
to the derivative of the Faraday tensor for which suitable field and evolution 
equations can be obtained — this idea is borrowed from the analysis of the 
conformal Einstein-Maxwell system of |B] . The system of equations considered 
are similar in spirit to those frame formalism equations of [7] . 

The evolution equations for the Einstein-Maxwell-Euler system were first 
considered from the perspective from the Cauchy problem in [5| — see also 
[3] . This analysis makes use of a traditional formulation of the Einstein field 
equations in which the gravitational field is described by the metric tensor, 
and uses wave coordinates to obtain non-linear wave equations for the metric 
tensor. In this approach, the resulting evolution system is of mixed order. As a 
consequence, the so-called Leray theory needs to be used — see As pointed 
out in that reference, a first order formulation of the problem does not nec- 
essarily improve the analytic results that can be obtained with a mixed order 
one. However, first order formulations may be preferred in numerical compu- 
tations and in the analysis of the motion of isolated bodies — notice, however, 
the local existence result for isolated bodies consisting of dust (charged and 
uncharged) given in [TO] which is based on a mixed order formulation. In any 
case, to the best of the author's knowledge, there is no frame formulation of 
the Einstein field equation with is used in current numerical simulations of 
the Einstein field equations. From our particular point of view, the value of 
first order formulations of the Einstein-Euler system like the one being used 
here lies in their amenability to a detailed analysis of their linear and non- 
linear stability — see pnil2] . Alternative discussions of the problem of the 
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well-posedness of the evolution equations of the Einstein-Maxwell-Euler sys- 
tem and of magnetohydrodynamics can be found in, for example, |13U14l[T51 
millT], and [H]. 

There are several formulations of the evolution equations of the Einstein- 
Maxwell-Euler system geared towards its use in numerical computations. With- 
out the aim of being exhaustive, we notice in particular [TOll^l^ — a good 
point of entry to the extensive literature in this respect is the review [35]. A 
peculiarity of the systems given in }191I201[^ is their explicit use of the con- 
ditions of ideal magnetohydrodynamics to construct their system. Hence, they 
may require some modifications if one is to use them to model more general 
types of electromagnetic fields. The question of the hyperbolicity of these sys- 
tems is by no means direct and has to be analysed with certain amount of 
care. In particular the reference [12] makes use of an ADM formulation for the 
Einstein part of the system which is known to have problems concerning hy- 
perbolicity — see e.g. [4]. This problem has been addressed in j20l[2T] by using 
a so-called BSSN formulation for the geometric part — see e.g. [23] for further 
details. The BSSN equations are first order in time and second order in the 
spatial coordinates and are known to satisfy, in the vacuum case, a certain no- 
tion of hyperbolicity (strong hyperbolicity) under certain additional conditions 
— see e.g. [23]. Local existence and uniqueness of solutions to strong hyperbolic 
systems follows under some further assumptions — see e.g. |25| . However, be- 
cause of the mixed order nature of the system, it is not directly obvious that 
this hyperbolicity property is preserved if one extends the vacuum Einstein 
system to the Einstein-Maxwell-Euler one despite the lower order coupling of 
gravity to the charged fluid. To the best of our knowledge this issue has not 
been addressed in the literature. A discussion on the caveats of this type of 
reasoning can be found, for example, in Section 8.5 of |25| . 

It should be pointed out that the evolution equations discussed in the 
present article are only suited to the discussion of unbounded self-gravitating 
charged fluids. In order to accommodate an initial boundary value problem, 
further structure needs to be verified in the equations. For example, it is not 
known whether our evolution equations are compatible with the class of max- 
imally dissipative boundary conditions used in, say |26| . to show the well- 
posedness of the initial boundary problem for the vacuum Einstein field equa- 
tions. In any case, it is interesting to point out that the analysis of |26| makes 
use of a frame formalism similar to the one used in the present article. 

Symmetric hyperbolicity, despite being a technical mathematical notion 
has a deep physical significance. Moreover, the question of the well-posedness 
of the evolution equations of the Einstein-Maxwell-Euler system is a problem 
that touches upon many aspects of current theoretical and numerical analysis 
of physical phenomena. Particular examples of these systems are given by by 
accretion disks around compact objects. Due to the complexity of the under- 
lying equations, the analysis of the stability and the shape of disks and flows 
is often addressed by numerical methods. Further applications of charged per- 
fect fluids in General Relativity can be found in, for example, [^[T71[?7ll28j . 
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Numerical simulations of magnetohydrodynamical systems with astrophysical 
applications in view have been a topic of interest in recent years — see for 
example [29.30,31,32,33]. 

The present article is structured as follows: The tetrad formalism used 
in this article is briefly reviewed in Sec. [21 In Sec. [3] we write and discuss 
the relativistic equations describing a charged perfect fluid. General remarks 
concerning the reduction procedure to obtain suitable evolution equations are 
given in SecU] the resulting evolution equations are discussed in Sec. 14.21 and 
the subsequent sections. The auxiliary field describing the covariant derivative 
of the Faraday tensor is presented in Sec. 14.31 A summary of the evolution 
equations is given in Sec. 14.61 The case of an infinitely conductive fluid (MHD) 
is briefly explored in Sec. 14.71 Some concluding remarks are given in Sec. [5] 
Finally, in Appcndix|X]we briefly discuss certain issues concerning homentropic 
flows and their relation to the equation of state of the perfect fluid. 

2 Tetrad formalism 

As mentioned in the introduction, our discussion of the Einstcin-Maxwell- 
Eulcr system will be based on a frame formalism. The purpose of the present 
section is, mainly, to flx the conventions to be used in the rest of the article. 

Let {Ai,gap) denote a spacetime. Our discussion will consider the two 
possible conventions for the signature of the metric tensor ga/s- We introduce 
frame fields {ea}a=o,...,3 on the spacetime satisfying gab = g[e.a,ei,) = 
r]ab = e diag(l, -1, -1, -1) and e = ±1 —see [MlIT] . We denote by (w") the 
corresponding dual basis (cobasisjl. The frame fields Cq and the cobasis w° 
are expressed in terms of a local coordinate basis as 

ea = ea"a„, Lu''=Lj\dx". (1) 

Thus, 

uj\eb" = S\, uj\ea^ = Sj. (2) 

The metric tensor and its inverse can be written as can be written in terms of 
and Ca as 

g^p = Vabt^^L^'p, .9"^ = TJ^^'ea'^eb^. (3) 

The coefficients of the commutator of the elements of the tetrad, D^^, are 
defined by 

eb] = DJe, = (ca^d^Cb^ - et^^^e/) Dp. (4) 

The connection components (Ricci rotation coefficients), -T^" > of the tetrad Cq 
arc defined by the relations 

VaCb = r^^e,, Vau'' = -ra^Lj'. (5) 

^ Latin letters a are the tetrad indices while Greek letters a denotes the tensorial char- 
acter of each tensor (spacetime indices). Latin letters {i, j, k, q,p), run in {1,2,3}, denoting 
quantities in the tetrad space. 
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In particular, since (rytic) — 0, one has the symmetry Fa^bc) — 0. The Rie- 
mann tensor^ ^"^bcd- satisfying the Ricci identity 

R\cd'^a = VcVd^a - VdVc^a - V[c,rf]Wh, (6) 

is given in terms of the connection coefficients i^j^" by 

In the sequel we will make use of the following decomposition of the Riemann 
tensor 

Rabcd = Cabcd + {ga[cSd]b - 9b[cSd]a} , (8) 

in terms of the Schouten tensor Sab 

Sab = Rab — ^Rgab, (9) 

and the Weyl tensor 

Cabcd = Rabcd + {ga[cRd]b - gb[cRd]a) + 3 9a[c9d]b- (10) 

In the previous expressions Rab = R'^acb denotes the Ricci tensor of gap , and 
R = g'^^'Rab its Ricci scalar. The components of the curvature tensor satisfy 
the Bianchi identity 

R\cd;e] = 0- (11) 

Now, from the contracted Bianchi identities for the Einstein tensor, Gab = 
Rab - \gabR-, satisfying 

V'Gab - 0, (12) 

we infer that 

Fabc = ^dFibc = 0, (13) 
where Fabcd denotes the so-called Friedrich tensor 

Fabcd = Cabcd - ga[cSd]b- (14) 

Taking the Hodge dual of Eq. (|T3)) with respect to the index pair (c, d) , we 
obtain an equation similar to for the tensor Fabcd defined by 

Fabcd = C'abcd + 2^pb^^acd^ (1^) 



where tabcd denotes the components of the completely antisymmetric Levi 

f 

cd ' 



Civita tensor with respect to the frame Ca and C*^^^ = \Cahef^^^ 
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3 Ideal magnetohydrodynamics (MHD) in curved spacetime 

In this section we introduce tlic basic equations describing a relativistic charged 
perfect fluid coupled gravity — the Einstein-Maxwell-Euler system. Tlie nota- 
tion and conventions used are based on those of [3] and p71l28j . 

We start by considering the following Einstein equations 

Gab = KTab, (16) 

with an energy-momentum tensor, Tab, which can be split as Tab = T^b + 
where 

Tib = (P + PWaUb - e pgab (17) 

is the energy-momentum tensor for an ideal fluid, while p and p are, respec- 
tively, the total energy density and pressure as measured by an observer moving 
with the fluid. The time-like vector field (flow vector) t/" denotes the normal- 
ized future directed 4- velocity of the fiuid. It satisfies 

U'^Ua = e. (18) 

Associated to U"^ we introduce the projection tensor 

hab = gab - f^UaUb, (19) 

projecting onto the three dimensional subspacc orthogonal to C/". Indices are 
raised and lowered with gab- Thus, one has that /i°f, = S°-b — (U°^Ub, h^ah"c ~ 
/i^, and h"-bUa = 0. 



3.1 The electromagnetic energy-momentum tensor 

The tensor T^J™ denotes the energy momentum tensor of an electromagnetic 
field: 

= (PacF, ^ - ^F,,F'--'gab^ , (20) 

where Fab is the electromagnetic field (Faraday) tensor. The latter can be split 
in its electric part, Ea = FabU^ , and its magnetic part, B"^ = ^e'^^'^'^UbFcd, with 
respect to the flow. Indeed, one has that 

Fab = ei2E^aUb] - eabcdB^W). (21) 

Alternatively, the one can write 

Fab — e {Fab — *Bab) , (22) 

where 

Fab = ^E^aUb], * Eab ^ ^abcdE'^U^ , Bab ^ '2B[aUb], * Bab ^ ^abcdB^V^ . 

(23) 

One can readily verify the properties 
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BabU'' = eBa, EabU'-^eEa, * BabU'' ^ * EabU'' = 0, (24) 

where, as before, * denotes the Hodge dual operator. In particular, one has 
that,**S„6 = -Bab, so that *Fab = e {*Eab - **Bab) = e {*Eab + Bab) , with 

Ba = *FabU'', *Fab ~ -^abcdF'^'^ , ** Fab = —Fab- (25) 

From the decomposition into electric and magnetic parts, one can readily verify 
that the Faraday tensor has six independent components as these vectors are 
both spatial — that is, EalJ"' = BalJ"' = 0. Using the decomposition into 
electric and magnetic parts, the electromagnetic energy-momentum tensor, 
Eq. (PO)) . can be written as 

T„7 = -'-UaUb{E^ + B^) + + B^) + Pab - 2eg^aUb), (26) 

where we have written E-^ = EaE°- and B"^ = B'^Ba, and Pab denotes the 
symmetric, trace-free tensor given by 

Pab = P^ab) = ^{E^+ B^) - {EaEb + BaBb). (27) 

Furthermore, 

Ga = eauvdE^B-"U'\ (28) 

denotes the Poynting vector. Alternatively, one can rewrite Eq. p6p in the 
form 

T^aT ^^{E^ + B^) - {EaEb + BaBb) - 2eg^aUb) - eUaUbiE^ + B^). (29) 

3.2 The Maxwell equations 

The Maxwell equations are given by 

V^aFbc] = 0, V'Fab = eJb. (30) 
Using Eq. one can rewrite them as 

V[a{Ebc]-*Bbc]) ^0, V'[Eab-*Bab)^Jb- (31) 

The homogeneous Maxwell equation can be rewritten in terms of the Hodge 
dual of Fab as 

V'"Fab = eV''{*Eab + Bab)=0. (32) 

Eq. d^T]) can be further expanded to yield 

[V^'Eb] Ua-Ea+ V''eabcdB'U^ + 2E[bV''Ua] = Ja, (33) 
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where E = V^V aEb stands for the covariant derivative of Ea along the flow. 
For Eq. p2p a similar computation gives 

V {eabcdE'U'' + 2B[,[/,]) = 0. (34) 

Projecting Eq. ([55)1 along the directions parallel and orthogonal to the flow 
vector one obtains the two equations 

eV'^Ea + Ebil'' - eabcdU"B''^'U^ = U^'Jt, (35) 

+h\,EaV''Ub - E,V''Ua = h\Jb, (36) 

Projecting Eq. along the directions parallel and orthogonal to U°' one 
obtains 

eabcdU^'E'^V'^U'' + eVBa - [/''Bb = 0, (37) 

tabcdh^fViE-'U'') + Bahl^fV^Ub - hI'fBb - BfVUa = 0, (38) 

where it has been used that V a{E^Ub) = EbVaU'' + UbVaE'' = 0. 

The electromagnetic current vector J° will be split with respect to the flow 
vector as 

J- = p^U''+f, (39) 

where pc denotes the charge density and is the orthogonally projected con- 
duction current. Using the decomposition Eq. p9p in Eas. psp and p6p. one 
obtains the following propagation equations for the electric and magnetic parts 



of the Faraday tensor: 

= -2E''hf[aVb]U'' - eabcdh'fV'^iB'^U'') - jf, (40) 

B^f^ = -2S''/iy[,Vfc]C/'' + eabcdh''fy''iE'U''), (41) 

and from Eas. ([57)) and ([55]) the constraint equations 

eD'^Ba = -eabcdU''E^V''U'', (42) 

eD'^Ea = eabcdVE'^V'-V' + epc, (43) 



where ihi^a) = hj'wb denotes the directional covariant derivative along the 
flow vector {Fermi derivative) and DaWb = ha^hb"\J uWv is the orthogonally 
projected covariant derivative of a vector. We note that Cabcdh f^V"^ {X'^U'^) = 
-curlX/ + eeafcdU^X^if^, where curlX/ = CfacdU'^VB". 

In the sequel, it will be necessary to specify the form of the conduction 
current, j". Consistent with Ohm's law, we assume a linear relation between 
the conduction current and the electric field. More precisely, we set 

f = <y'"'Eb, (44) 
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where cr"'' denotes the conductivity of the fluid (plasma) . We will restrict our 
attention to isotropic fluids for which cr"'' = ag"^^, so that Eq. (j39p becomes 

J" ^ PcU" + (7E'', (45) 

with (T the electrical conductivity coefficient. Ideal MHD is characterised by 
the condition, cr — ?> oo (ideal conductive plasma). It follows from Eq. (|45|) that 
for this to be the case one requires the constraint Ea = 0- 



3.3 The fluid equations 

Form the conservation of the energy-momentum tensor 

yTab = 0, (46) 

and using Eqs. ((T7)) and ([^0]) and the Maxwell equation Eq. ([50)) one obtains 

V"TX = -e(V'^F,,)f^,^ (47) 
V°Tf, = UbUaV'^ip + p) + ip + p) [C/b(V"C/J + C/aV"C/b] - eVfcp. (48) 

Therefore, Eq. reads: 

C/6C/aV"(p+p)+(p+p) [[/^(V'^f/,) + t/aV^t/bl-eVbp-e (V"F„,) ^^^^ - 0. (49) 

In what follows, we will consider the projections of along the directions 
parallel and orthogonal to the flow lines of the fluid. Contracting Eq. with 
ell^ we obtain the conservation equation 

t/aVV +{p + p)V''Ua - U'^F^^iVFac) = 0, (50) 

where the constraint pSD has been used. For an ideally conducting fluid, where 
Ea = FabU^ = 0, the last term of Eq. ([50)) is identically zero and the electro- 
magnetic field does not have a direct effect on the conservation equation along 
the flow lines. Contracting ([^^ with the projector ft,'"^ one obtains the Euler 
equation 

{p + pWVaU^ - eh'-VbP - e{V''Fad)F,''h''- = 0. (51) 
The later equation can also be written as 

(p + ppaV'Ub + {UbV'Vd - eVb)p - €{V''Fad) {F,^ - eF'''U,Ub) - 0. (52) 

In the ideal MHD case the last term of Eq. ([52)) is identically zero so that 
Eq. ([5T)) reduces to 

{p + pWVaU'' - eh'"'\/bP - e{V''Fad)F'"' = 0. (53) 
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3.4 Thermodynamic considerations 

Here we consider a one species particle fluid (simple fluid). We denote by 
n, s, T the particle number density^ the entropy per particle and the absolute 
temperature^ as measured by comoving observers. We also introduce the volume 
per particle, v, and the energy per particle, e, via 

v = -, e = ^. (54) 
n n 

In terms of these variables the first law of Thermodynamics, de = —pdv + Tds, 
takes the form 

dp = ^-i-^drt + nTds. (55) 
n 

Assuming an equation of state of the form p ~ f{n, s) > 0, one obtains from 
Eq. dSSl) that 

p{n,s) =n(^^^ ~p{n,s), ^("' = ^ (^|^) ■ (56) 
Assuming that dp/ dp > we define the speed of sound, Vg = Vsin, s), by 

-?-(|^) ~|^>0. (57) 
\op J p + pon 



Since we are not considering particle annihilation or creation we consider 
the equation of conservation of particle number: 

V^VaU + nVaU" = 0. (58) 

Combining this equation with Eqs. (|50|) and (|55p we obtain 

U'''7aS=^U'F,-V''Fac. (59) 



In the case of an infinitely conducting plasma, where the last term of Eq. (j55p 
vanishes, Eq. ((59)) describes an adiabatic flow — that is, If^VaS = 0, so that 
the entropy per particle is conserved along the flow lines. A particular case of 
interest is when s is a constant of both space and time. In this case the equation 
of state can be given in the formp = p(p)- A further discussion on homentropic 
flows and barotropic equations of state is provided in Appendix |21 
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4 The hyperbolic reduction procedure 

4.1 General considerations 

Following [3], we introduce the notation 

iV" = <5o^ iV EE N-ea = eo. (60) 

For a given tensor any contraction with N will be denoted by replacing the 
corresponding index by N , and the projection with respect to hg'^ will be 
indicated by a prime. Accordingly, for a tensor Tabc one writes 

Tl^, = T^pgha'^NPhb". (61) 

Introducing the notation Eatc ~ ^Nabc^ where £0123 = 1, we obtain the decom- 
position 

Cafccrf 2e {Nlaeb]cd - <^ab[cNd]) ■ (62) 

Given a spatial tensor Ta-^...ap = T'ai - ap define the spatial covariant deriva- 
tive 

VaTa,...a, = ha''ha,'" ■ ■ ■ K^^'-V bTb^.-.b^. (63) 

One has. in particular, that T^ahbc = T^a^bcd = 0. 
For convenience, we introduce the tensors 

a'' = 7V''Vfc7V^ Xab = KV,Nb, X = h'"'xab. (64) 

In terms of these we can write 

VaTV' = eiV.a^ + Xa', a" = hT^o. Xab ^ -eha^h'^T^d- (65) 

Key for our subsequent discussion will be the decomposition of the Weyl 
tensor on terms of its electric, Eab = C^aNb^ ^^'^ magnetic parts, Bab = 
^* NaNb^ with respect to Na given by 

Cabcd = 2e(/t[ci?d]a - K{cEd]b) - 2(iV[ci3d]pe^ab + N[aBb]p(-^'cd) (66) 

where lab = hab — ^NaNb- From the Bianchi identity Eq. ([T^ for the tensor 
Fabcd we obtain the decomposition 

Fabc = Na {F'^bNNc ~ i^^cW^b) ^ '2eF^N[bNc] + e^a^^W + ^ahc, (67) 

from where it follows that 

'^F'aNb ~ ^NF'j^aNb + ^^'^F^aNb ^ ^'^ {F'Nacb + F'^aNb) + ^'^a-F'wArWb 

-^X" ■P'cadb - XaF'NcNb ~ Xb^F^aNc + X^aFcNNb + XF'^aM^) 
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4.2 Detailed analysis of the evolution equations 

In the sequel we make the following specific choice of tinielike frame vector; 

N^eo = U, Na = Ua. 

Using the decomposition given by Eq. (|14p we compute the following com- 
ponents of the tensors Fated and Fabcd- 



^NNNa 



^NaNb 



^aNNb 



^Nabc 



0, F' 



NNNa 



J^ab, NaNb ~ ^ab + "g~ 2" °' 



^aNNb = ^Bab + -Z {'^p^U'") £aub^" ^ 



-Eab + l^^ab 



(69) 
(70) 
(71) 
(72) 

(73) 



^Eap^l^ —SaubcU^ + {Tp™h" a) e^ufccC^", F'Nabc — ^fBape^bil^) 



^aNbc 



^abNc 



^abNc 



pi 

abed 



pi 

abed 



en /2 



Kme = -eEape^be + y ( 3P + P ) Up^'abe + ^ (T^^C/") d. fh\h\hf e. (75) 



Bape^be - % [hab {Tf^V'^h^) ~ Ke {T]'^U^hfb)] 



^ JJV I l^rpcmiu p TJV 



2e 



(76) 
(77) 
(78) 

-Vabe'^ed - ^ep«e//l^/l"a/l^c/l^rf + | {T;^'h\) e\ef h\h^ eK^Q) 

(80) 



— —^BcpC^ab + -zhacT°^U^h^'b, 



h[eEd]a — lalcBd]b 



hc[ahb]d ^ i^T°^K" hhr\„h, 



uv b'*c[a"d] 



Moreover, Eq. © becomes, 



^ab — Oab + ^ ab 



(81) 



where 



'ah — 



ry-if edrjii 

^ab ^ '^yabg ^ ed 



(82) 



and 



'-'ab = '^^ah 



-en [ FaeF'^'gdb - ^^^cdf-^Vb 



(83) 



Using Eq. (j8ip and Ec^. ((66|) in Eq. ([8]), the Riemann curvature tensor Rabcd, 
can be decomposed in the form 



Rabcd — R™bed + ^abed + ^a' 



bed^ 



(84) 
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where 



and 



(.9a[c.9d]6 ~ 9b[c9d]a) ~ '^(-U[bga]lcUd] 



+2e [U[agb][cQd] - G[a9b][cUd]] 



(85) 



K^bcd = '^^{h[cEd]a - L[cEd]b) - HU[cBd]p<^''ab + ^[aBb]p^^ cd) ^ 
B-abcd = '^'^'^P {^b[chd]a " h[chd]b) + K-P {ha[cUd\Ub - hb[cUd]Ua) 



(86) 
(87) 



Finally, from the decomposition Eq. 



^^(a|Ar|h); 



we identify the tensors 
1 . 



■^(a|JV|6) - -^^abh^^^Fl^ 



Nv 



As it will be shown in the next subsections, the evolution equations for the 
electric and magnetic parts of the Weyl tensor will be extracted from these ten- 
sors. The essential difficulty in our analysis resides in a satisfactory treatment 
of these equations. 

4-2.1 Evolution equation for Bab 



From Eqs. ([BS]) and (|13p we obtain the following evolution equation for the 
magnetic part of the Weyl tensor 



= e^(a|Ar|b) — C-uBab - D dE c(a^b)'^'' + '^^O'C^'^'^ (aEb)d - X'^{aB: 



b)c 



-2x(a^Sfc), + xBab - eXcdBpqeP\y\^ + eF^l 



where 

^-^(alAflft) 



{a\N\ 



(89) 

In vacuum, Eq. (|88p together with its electric counterpart provide hyperbolic 
evolution equations for the electric and magnetic parts of the Weyl tensor 
independently of the gauge choice. This is not automatically the case in the 
presence of matter. Using Eq. ([55)) in Eq. we find that the electromagnetic 
contribution to the evolution equation of Bab is given by 



P em 

^-^ia\N\b) 



~2~ 



cv{b 



"^A uef 



h^bh" a)F.cFj,^ - \hf ^i,ha)pFq.F'^'^ - U^eP,.aibX\)EdFp^) 



Notice that this last expression contains derivatives of the Faraday tensor 
which cannot be replaced by means of the Maxwell equations. These deriva- 
tives enter into the principal part of the evolution equations and destroy the 
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hyperbolicity of the evolution equations for the magnetic part of the Weyl 
tensor. In order to deal with this difficulty and additional variable, represent- 
ing the derivative of the Faraday tensor needs to be introduced. This will be 
discussed in subsection [ 



4-2.2 Evolution equation for Eab 

The evolution equation for the electric part of the Weyl tensor can be obtained 
from Eq. (j68[) as follows: 



= e 



^UN\b) 2 



— CljEab + DcBd{a^b)'^'^ — '^^O'cf^'^ {aBb)d " "^Xia^ Eb)c — '^X'^iaEb) 



cd 



+ habX Eed + 2xEab + ^{p + p) 



X(ab) - -^Xhab 



TTi'cm 

^(a\N\b) - 2 



T:habh^''F:S'M) 



where 



p ern " J uuv 771 cm 

^(a\N\b) ~ -^ii-abn' ^uNv 



\eKCu {Trh\h\) + y 



rpCniT rUTjV 



hnb 



X{ab) 



2X(a h b)- X(a rib) 7, — 



-X 



(92) 



Finally, using Eq. (j83p we obtain the explicit expression 



P em 
ia\N\b) 



habh''" 



- uNv 



kCt 



IT- (an b)^qf-fc 



habV' 



{E^Fy,) 



-KeE^Fyc [ h^iaO-b) ~ -^a" 



-F ,iF " 



2X(a hb) - X {ah) ~ 2^ 



-1x{ab)FqpF''P - y£;\(a,,) " ^xhabF^pF^^ + ^eKxhabE^ 



(93) 



As in the previous subsection we observe the presence of derivatives of the 
Faraday tensor which need to be dealt with by the introduction of a new field 
if one is to preserve the hyperbolicity of the equations. 



4.3 An auxiliary field 

As discussed previously, derivatives of the Faraday tensor appear in the evo- 
lution equations implied by the Bianchi identity. This feature of non- vacuum 
systems destroys, in principle, the hyperbolicity of the evolution equations. 
In order to get around with this difficulty, it is necessary to introduce the 
covariant derivative of the Faraday tensor as a field variable. More precisely, 
let 

^abc = yaFbc- (94) 
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As a consequence of the definition and tlie Maxwell equations the tensor t/^abc 
has the symmetries t/jabc = 4'a[bc] , V'[abc] = 0. In order to construct a suitable 
equation for ipabc "we consider the commutator of the covariant derivative V in 
the form 

VaVfcF,rf - VbVaFcd = ~2F,[,R''d]ab- (95) 
Making use of the definition of tpabc one obtains 

Va^Abcd - Vb?Aacd = -2F^[cR''d]ab, (96) 

where in the right hand side we use the expression for Rabcd involving the its 
decomposition in terms of the Weyl tensor and the matter fields. The tensorial 
expression (|96[) provides the required field equations. 



4.4 Propagation equations for the auxiliary field 

In order to deduce from them suitable evolution equations one proceeds by 
analogy to the case of the Faraday tensor Fab- Proceeding as in the case of 
Eqs. pop and ([5^ we obtain the following two equations for ipabc 

V^'i^adb = 2F'-^''Rd^aeb - eVa Jd, (97) 
V-^,,, = er'F^uR'^a,, (98) 

where Rabcd is given by Eq. (|84p . In the subsequent discussion is convenient 
to introduce here the following notation 

^l;ayz='l'yz='F[yz]- (99) 

That is, the first index of ipabc tensor is dropped. Because of the ipabc is an- 
tisymmetric with respect to the last pair of indices, it is natural to introduce 
its electric and magnetic part respect to Ua via 

£d = 'FdN = i'adnU", = ie^-'C/^tf'.t, (100) 

respectively, so that are £d and Bd are spacelike vectors: 

SaU'' = 6aC/° = 0. (101) 

Accordingly, we write the tensor •f' as 

1'ab^e[£ab-*Bab]- (102) 

Furthermore, one has that 

*'^a6 = 6 [*£ab + Bab] , (103) 

where 

£ab = 2£[aUb], *£ab = ^abcd£'^U'^ , (104) 

Bab^2BlaUb], *Bab = eabcdB'U''. (105) 
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Continuing the analogy with the Maxwell equations and the Faraday tensor 
Fab, it follows that Eqs. and Eq. read 

£d - (V^fb) Ud - V^'edbcvB'U" - 2f[fcV^C/d] = e5d, (106) 

V {eabcdS^W + 2B[aUb]) = eVb, (107) 

where 

Sd = '2F'^^''Rd]feb ~ ^'^ fJd, = —'^f^'^FeuR'^dab- (108) 

Projecting Eq. p06p along the directions parallel and orthogonal to the 
4- velocity U^, one obtains the following set of two equations: 

Sbil'' + eV-'Ea - eabcdU''B''V''U'' = -eU'^Sb, (109) 

+V"'{eaecdB'U-') = -eh\Sb. (110) 

A similar procedure applied to Eq. (jl07p gives 

EabcdU^'S^V'^U^ + eV^'Ba + B'Ub - eU'Vb, (111) 
ea&cd/i''/V''(rC/'^) + Bah''fV''Ub - h^fBb - BfV'^Ua = e/i/Vb. (112) 

The propagation equations p09[) - (jll0p can be written as 

- -2£''hf[aVb]U' - eabcdh'fV^B'^U'') + eh^fSd, (113) 
B^f-^ = -2B''hf^aVb]U^ + eabcdh^fV''{£''U'') - eh^' fVd- (114) 

while the constraint equations (|llip - (jll2p assume the form 

eD-'Ba = -eabcdU'^E'VU-' + eU-'Vd, (115) 
eD"-£a = tabcdU'^B^'V^U'' - eV'Sd- (116) 

In order to verify that Eqs. (jll3p - (jll4|) give appropriate hyperbolic evolu- 
tion equations, we need the explicit expression for the electromagnetic current 
vector J" given by Ohm's law — see Eq. ps)) . A computation then gives that 

Vbh^v = —^f^^FeuR'^dab^v- (H^) 

Sdh'', = F^'-Rdfebh^ + Rfe (n - ^E^Uv) ~ e {V f{oE,) + eU^aEdV + PaV fU^) . 

In the last equation it is assumed that the derivatives of the components of the 
Faraday tensor appearing in the right hand side of (|118p are rewritten in terms 
of the tensor V'abc- Furthermore, Rabcd is expressed in terms of its irreducible 
components. 
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Finally, using the kinematic quantities defined in Eq. (j64p we finally rewrite 
Eqs. (|ll;/flll4l) as: 

^(/> = [xha - hfb {eUaa'' + x'a)] + cuTlBf - eeafcdB'U'^a'' + eh"" fSd, 

(118) 

B^f) = -2B^ [xhfa - hfb {eUao!' + x\)] - curlf/ + eeafcdS'U^a" - eh^fVd, 

(119) 

where Vbh''y and Sdh'^v are given by (|117p - (|118p . These equations provide the 
required evolution equations for the components of the auxiliary field ipabc- The 
hyperbolicity follows by comparison with the Maxwell evolution equations, 
Eqs. (Unil-dlll). 

4.5 Propagation equations for frame coefficients and connection coefficients 

We now consider the evolution equations for the tetrad coefficients, ea^, and 
the Ricci rotation coefficients, /l,°c- 

The gauge choice N — eo — U, Na — Ua, = eg" ~ Sq"^ together with 
Eq. (HI) readily give 

5*e/ = (V,~r/o)e/ + ro%^ (120) 

In order to obtain evolution equations for the Ricci rotation coefficients, 
one imposes the gauge condition Fo^i = arising from the Fermi propagation 
of the spatial frame vectors — see [3] . It can be verified that as a consequence 
of the metric compatibility condition -Tajfjc) = one only needs to consider 
evolution equations for the component F^^j, Fq"^ and 1^^°. 

A propagation equation for F^^^ is readily found from Eq. ([7]). One has 
that 

The propagation equations for Fq'^^ and F^^^^ require a more elaborated 
analysis involving the fluid equations. From Eq. together with the ther- 
modynamic relations (|5^ we obtain 

p = -l^y^Uaip + p) - U'^F^^V'^Faa. (122) 

Furthermore, using Eq. ((50)) and Eq. (jl22[) we have 

eVdP ={p + p){Ud - UdvlV'Ua) - eV'^FacF/, (123) 
so that the condition Vj^V^jp = with Eq. p23p implies 

= Ve [ip + p){Ud - UdvlV'^Ua)) -Vd ({p+p){Ue " U^vlV^Ua)) -2e (^^[j=Ve] + V[eFj;) VF, 

(124) 
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This equation can be conveniently split as the sum of the three terms: 

Jed = Aed + Med + J^ed = 0, (125) 

where 

Aed = 2ip+p)V[eiUd] - Ud^ylV'^Ua), (126) 

Med = 1{U^d - f/[d'^fV''[/|,|)Ve](p+p), (127) 

:Fed = -2e (Fj^V,] + V[,F^f ) V^Fae- (128) 

Now, if one regards the velocity of sound Vg as function of n and s, then from 
the (0, i) component of Eq. (|125p and Eq. (|123p we obtain: 

f S ) ^/o f^o". + ^{v'Fde)F, ^) - ^r,\ Vo. + ^v..r/o 
^; \"P J s \ p + p J p + p p + p 

- "^^'^"^ " ^V[o(^^,]^V'^F,,), (129) 

where 

la = '^^^F^'V'Fdc, (130) 

Similarly, from the (j, fc) component of Eq. (|125p we find that 

OR 2 2 

+ ^^[^^^0^1 + ^^o%7,] + ^V[,(F/V'^^^..) = 0. (132) 

Finally, making use of Eq. (jl2ip in Eq. (jl29p one obtains 

+ {r,\ + ^{v''Fde)FA -^r,\Vos + ^v.sr/o - ^(r„",7o - r^^i.) 
\ p + p J p + p p + p p + p 

--|-V[o(F,]^v'^F,,) + (i?^o,, + r,^o{r, ° - r^\) - r/or/, + r/^r/o) , (i33) 

while using 

e^ro",- + i?°o;, + ro° ro"^, - F, ° rj,'o - r^.^^ (134) 
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from Eq. (Q in Eq. (fT^ gives 
2 2 

(135) 

In all the previous it is understood that terms of the form V^Fdc are to be 
replaced using the relation 

V'^Fdc^i'c = PcJJc + (tE,, (136) 

where it is assumed that the eharge density, pc, is constant multiple of the 
fluid density, p = gpc ■ 

Remark. Important for the hyperbolicity of the evolution equations is that 
the term 

-^V[.(nfV^^^..), 

appearing in Eqs. (|129p and Eq. (|132|) . can be rewritten in terms of (-T^'d' •^a, Sa) 
if one uses Eqs. (|TUgi) and ([TM)) . 



4.6 Summary of the analysis 

In this subsection we summarise our analysis of the evolution equations. Evo- 
lution equations for the independent components of the vector variable 

V= (^ei'^'',rQ°i,m,Eab,Bab,Ea,Ba,£a,t3a,n,S,Sa^ , (137) 

have been constructed. The components of il-'fc not included in this list are 
determined by means of gauge conditions and symmetries. By construction, 
the electric and magnetic parts of the Weyl tensor are tracefree. This symmetry 
is disregarded in these considerations and is recovered by imposing it on the 
initial data. It can be shown that if these tensors are initially tracefree, then 
they will be also tracefree for all later times — see e.g. [5]. 

The evolution equations for the independent components of the unknowns 
in (|137p and the underlying assumptions arc given as follows: 

(i) The propagation equation for the tetrad coefficients, e^'' , is given by Eq. (|120|) 
by virtue of the Lagrange condition [/° = eg ~ Sq . Eq. (|120p has the same 
principal part than the corresponding equation in the uncharged case as 
described in [3] . It gives rise to a symmetric hyperbolic subsystem of equa- 
tions. 
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(ii) The connection coefficients, 1],° and I^j"^ are given, respectively, by Eqs. 
(|129p and Eq. (|132p . In addition, the gauge condition, Tq-'- = 0, is assumed. 
Eq. (|121[) takes care of Fj \ . As in the case of the equations frame coeffi- 
cients, Eqs. (jl29p . (jl32p and (|121[) have the same principal part as in the 
uncharged case analysed in [3]. These equations again form a symmetric 
hyperbolic subsystem of equations. 

(iii) The evolution equations for the electric part, -Bah, and the magnetic part, 
Bab, of the Weyl tensor arc given, respectively, by Eqs. (|91|) and ([55)) . As 
mentioned before, the traccfrecness of these tensors is not used to reduce 
the number of independent components. Thus one has 12 equations for 
equal number of components. Equations with a principal part of the form 
of Eqs. (|?T|) and ([55)1 are symmetric hyperbolic independently of the gauge 
choice — see e.g. [5]. 

(iv) The propagation equation for the electric part, Ea, and the magnetic part, 
i?b, of the Faraday tensor are given, respectively, by Eqs. (|40)) and (|4T|) . 
The isotropic Ohm law, Eq. (^5]) has been assume in the deduction of these 
equations. As in the case of the equations for the electric and magnetic 
parts of the Weyl tensor, the principal part of these equations is known to 
be hyperbolic independently of the gauge — again, see e.g. [2]. 

(v) The evolution equations for the electric part, f^, and magnetic part, Bf,, 
of the field V'a&c, corresponding to the covariant derivative of the Faraday 
tensor are given, respectively, by Eqs. (jllSp and (jllOp . These equations 
involve 24 equations for as many unknowns. Their structure is analogous 
to that of Eqs. (^0]) and (HU, except that they contain an extra free index. 
As a consequence, their principal part gives rise to a symmetric hyperbolic 
subsystem. 

(vi) The evolution equation for the particle number density n is given by 
Eq. (|58p . This equation is consistent with symmetric hyperbolicity as in the 
Lagrangian gauge all the derivatives of the flow vector U"' are replaced by 
connection coefficients, and thus, it contains no derivatives in the spatial 
directions. 

(vii) The evolution equation for the entropy per particle, s, is given by Eq. (j59p 
with the understanding that the covariant derivative of the Faraday tensor 
arising in the right-hand side is to be expressed in terms, either, of the 
auxiliary field ipabc or written in terms of the current vector Ja using the 
inhomogeneous Maxwell equation (pO|) . 

(viii) Finally the equation for Sa = VaS is obtained by covariant differentiating 
Eq. (jl30p . commuting covariant derivatives and taking into account having 
in mind Eqs. (|136p and, again, (|130p . 

Summarising, Eqs. (gOl), dUD, dUl), (HH]), (HSl), 1^,, 

(jl29l) , p32p and (|130|) provide the desired symmetric hyperbolic evolution sys- 
tem for the Einstein-charged perfect fluid system. This system can be written 
in the form 



A°5ov - A^djX ^ Bv, 



(138) 
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where A'^ — [x^ , v) arc matrix- valued function of the coordinates and the 
unknowns v, and A*^ is positive definite. The structure of the system (|138p 
ensures the existence of local solutions to the evolution equations. The analysis 
of whether these solutions give rise to a solution of the full Einstein-charged 
perfect fluid system requires the analysis of the propagation of the constraint 
equations. This is a computationally intensive argument which will be omitted 
here. 

4.7 Infinitely conductive plasma 

A particularly important subcase of our previous analysis is that of an ideal 
conductive plasma, defined by the condition Ea = — ideal magnetohydro- 
dynamics. The hyperbolic problem for this case can be naturally treated as a 
subcase of the general case in which the Faraday tensor has all the six inde- 
pendent components. In this case the vector variable v reduces to 

V = {ef , 1^0° . ^» Eab, Bab. Ba, £a, Ba,n, S, Sa}. (139) 

Notice that although Ea = 0, one nevertheless has a no n- vanishing electric 
part of the auxiliary tensor ipabc as can be seen from £i = "Bj^Fijlf^. 

In the infinitely conductive plasma case the system consisting of Eq. (j30[) , 
Eq. dSO]), Eq. and Eq. ([59]) reduces to 

t/aV>+(p + p)V''[/, =0, (140) 
{p + pWVaU' - eh'^-'VbP - e{V''Fad)F,''h'^ = 0, (141) 
U^VaS = 0, (142) 

Remark 1. As a consequence of Eq. (|142p the entropy per particle is constant 
along the flow lines of the perfect fluid. This fact is translated into Eq. p42p , 
propagation equation for s, that implies the propagation equation CjjSa = 0. 

Remark 2. Consistent with Ohm's law the source term of the Maxwell equa- 
tion is in this case simply given by Ja = PcUa- 

5 Conclusions 

In the present article we have revisited the issue of wellposedness initial value 
problem for the evolution equations of the Einstein-Maxwell-Euler system (a 
selfgravitating charged perfect fluid). Our analysis assumes a system consist- 
ing of a single relativistic perfect fluid obeying to a certain equation of state. 
The approach followed makes use of the well known 1-1-3 tetrad formalism by 
means of which the various tensorial quantities and equations are projected 
along the direction of the comoving observer and onto the orthogonal subspacc. 
Following , we require the timelikc vector of the orthonormal frame to fol- 
low the matter flow lines (Lagrangian gauge). Moreover, we assume the vector 
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fields tetrad to be Fermi transported in the direction of U, tliese conditions 
fix certain components of tfie connection. 

A key feature of our analysis was the introduction of a rank 3 tensor tpabc 
corresponding to the covariant derivative of the Faraday tensor Fab- The pur- 
pose of introducing this tensor was to ensure the symmetric hyperbolicity of 
the propagation equations for the components of the Weyl tensor. The sym- 
metries of the tensor ipatc readily allow a decomposition analogous to that of 
the Faraday tensor into electric and magnetic parts. 
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A On the equations of state 

A.l Homentropic fiows 

A fluid is said to be homeontropic is s is constant in space and time. In general, the equation 
of state is given by p = f{n,s). Now, if one has an homeontropic flow, the latter can 
be written as p = /(n) — there is no dependence on s as it is constant. Now, if / is a 
differentiable function of n and f'(n) ^ 0, then one can write n = /~^(p). It then follows 
that 

dp 

V = n- p, 

an 

= nf'{n) -p, 

= nf'{r\p))-P- 

Thus, p is of the form p = h(p) — that is, one obtains a barotropic equation of state. 



A. 2 Conditions for a barotropic equation of state 

In the previous section we have shown that the assumption of an homeontropic flow implies 
a barotropic equation of state. We now investigate more general situations for which one 
can have this type of equation of state. 
Assume that p = h{p) . Then, from 

p{n,s) = nf^\ -p{n,s) 



it follows that 



M/(n,s)) = n(^|^) -/(n,s). 
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The latter can be read as a (possibly non-linear) differential equation for /(n, s). It can be 
integrated to obtain 

where g{s) is an arbitrary function of s. This expression can be used to eliminate n from 
the discussion. In the case of a dust equation of state (p = 0), the last relation reduces to 

p/n = g{s). 
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